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The superconducting flux qubit has two quantum states 
with opposite magnetic flux. Environment of nuclear spins 
can find out the direction of the magnetic fiux after a deco- 
herence time To inversely proportional to the magnitude of 
the flux and the square root of the number of spins. When 
the Hamiltonian of the qubit drives fast coherent Rabi oscil- 
lations between the states with opposite flux, then flux direc- 
tion is flipped at a constant rate a; and the decoherence time 
T = ujtq is much longer than tq. However, on closer inspection 
decoherence actually takes place on two timescales. The long 
time r is a time of full decoherence but a part of quantum 
coherence is lost already after the short time tq. This fast 
partial decoherence biases coherent flux oscillations towards 
the initial flux direction and it can affect performance of the 
superconducting devices as qubits. 

PACS: 03.65.Yz, 03.67.Pp, 85.25.Dq 

Quantum computers can perform certain tasks expo- 
nentially faster than their classical counterparts [1] . The 
basic unit of the quantum computer is a quantum bit 
or qubit. Several physical implementations of the qubit 
were proposed including ion traps [2], nuclear magnetic 
resonance [3], quantum dots [4], single photon [5], and 
two flavors of superconducting Josephson junction qubits 
[6] . The solid state implementations arc in principle most 
easily scallable and most compatible with the existing 
classical electronics, but at the same time they are the 
most endangered by quantum decoherence. However, re- 
cent experiments [7,8] in both types of superconducting 
qubits show hundreds of coherent oscillations. This is still 
far from the minimal requirements of the fault tolerant 
quantum computation [9] but it makes superconducting 
qubits attractive candidates for scallable quantum tech- 
nology. In order to push their performance even further 
it is essential to understand their quantum decoherence 
better. 



I. SUPERCONDUCTING FLUX QUBIT 

The superconducting flux qubit is a superconducting 
loop broken by one or more Josephson junctions [6]. The 
junctions are designed so that the lowest two quantum 
eigenstates are states with opposite magnetic flux. The 
flux can range from a tiny 10"'^ . . . 10~^ fraction of the 
flux quantum [10] to 1/2 of the quantum [11]. The two 
states span a two-dimensional Hilbert space of the qubit. 
I assume convention that the state |1) (|0)) with posi- 
tive(negative) flux is the -M(— 1) eigenstate of Cz- Degen- 
eracy of these two states is removed by coherent tunnel- 



ing between the opposite flux states driven by the qubit 
Hamiltonian: 



(1) 



For /i = the eigenstates of Hq are the coherent superpo- 
sitions of opposite magnetic fluxes |±) = ^^^^^ ■ These 
Schrodinger cats were observed in two independent ex- 
periments [10,11]. For /i = the qubit Hamiltonian (1) 
drives coherent oscillations between the states with op- 
posite flux. For example, the initial state |1) evolves into 
|1) cos{\u)t) — «|0) sa\[^ujt) with an expectation value of 
the flux operator 



(2) 



oscillating with the Rabi frequency oj. Damped flux os- 
cillations were observed in a number of experiments [7,8] 
in both types of the superconducting qubit. 

The qubit Hamiltonian (1) is sufficient to perform any 
one-qubit operation. Two-qubit operations are possible 
thanks to inductive coupling between any two qubits [6] . 
Entangled eigenstates of two qubits were detected in Ref. 
[12] and coherent Rabi oscillations in a system of two in- 
ductively coupled qubits were observed in Ref. [8]. It 
is possible to perform single-qubit NOT and two-qubit 
CNOT operations in an adiabatic way [13] that, in par- 
ticular, precludes excitation beyond the truncated two- 
dimensional Hilbert space of the qubit. However, before 
it comes to quantum computation quantum decoherence 
has to be overcome first. 



II. QUBIT IN A SPIN BATH 

Various sources of decoherence for flux qubits have 
been discussed in Ref. [14]. At relatively high temper- 
atures the main source of decoherence are normal state 
quasiparticles. However, density of quasiparticles is ex- 
ponentially suppressed at temperatures much less than 
the critical temperature. Other decoherence mechanisms 
including electromagnetic radiation from the qubit or 
ohmic dissipation in the environment are respectively 
negligible or tractable. Nuclear spins are argued to play 
a minor role. Due to relaxation the spins randomly flip 
their polarization. Random spins are a source of random 
magnetic field which couples to the magnetic moment of 
the qubit and randomizes its quantum state. This pic- 
ture is further corroborated in Ref. [16] where the spins 
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are assumed to be mutually non-interacting but each of 
them is coupled to a bosonic environment. The dccoher- 
ence time is estimated in the range of miliseconds. The 
final Eq.(47) in Ref. [16] shows that in the limit of van- 
ishing nuclear spin relaxation rate (7^ in Ref. [16]) the 
decohcrence time tends to infinity. As expected, the ex- 
ternal magnetic noise vanishes for vanishing nuclear spin 
relaxation rate. 

However, apart from generating magnetic noise nuclear 
spins can be also silent witnesses of the quantum state 
of the qubit. In the limit of vanishing spin relaxation 
and for negligible spin-spin interaction each spin is sim- 
ply precessing in the magnetic field of the qubit. The 
direction of the field and the direction of the precession 
depend on the state of the qubit. This way the spins can 
learn the quantum state of the qubit. Once they know 
the state any quantum coherence between the states with 
opposite flux is lost. This elementary argument shows 
that decoherence does not vanish for vanishing spin re- 
laxation. 

In this paper I neglect nuclear spin relaxation. When 
the relaxation rate is longer than any other relevant 
timcscale, like the frequency of the Rabi oscillations or 
frequency of the spin precession in the magnetic field 
of a qubit, then this is a very reasonable assumption. 
Spin lattice relaxation is also neglected in the paper of 
Prokofev and Stamp [15]. In contrast to the very general 
formalism employed in Ref. [15], in this paper I apply 
most elementary methods to a simple model Hamilto- 
nian. This simple approach benefits with clear interpre- 
tation of the results and additional insights into dynam- 
ics of the decoherence process. In particular, the sim- 
ple model makes it very clear that even in the absence 
of any spin lattice relaxation the spin environment deco- 
heres quantum state of a single qubit. The decoherence is 
not an artifact of ensemble average over different static 
spin configurations or an ensemble of qubits, as some- 
times claimed in the literature, but a result of genuine 
entanglement between the qubit and the spins. 

The qubit interacting with spins is described by a 
Hamiltonian 

H = Hq + V + Hs . (3) 
Here V is interaction between the qubit and N spins 



N 



N 



(n) 



(4) 



n=l 



B{fn) is a magnetic field of the qubit in the state ]1) at 

the position r„ of the n-th spin. Direction of the mag- 
netic field is reversed in the state |0) as is accounted 
for by the operator (T^. ct^"^ is a vector of Pauli matrices 

^ci"^ , (Ty"^ , (7^"^^ in the Hilbert space of the n-th spin. A 
unitary transformation in the Hilbert space of each spin 
brings the interaction Hamiltonian to a more convenient 
form 



(5) 



Here i3„ = |i?(r„)| > is the strength of the qubit mag- 
netic field at the location of the n-th spin. 



III. EXACTLY SOLVABLE MODEL 

When the magnetic field from the qubit is stronger 
than magnetic fields from other spins, then spin-spin in- 
teraction can be neglected, Hs = and 



H 



(6) 



This exactly solvable spin-spin model was also considered 
in Ref. [17]. The Hamiltonian (6) can be easily diagonal- 
ized. The spin part of any eigenstate is 



\s) = \si) . . . |sjv) . 



(7) 



Here |s„) is an eigenstate of cri"' with an eigenvalue s„ e 
{-1-1,— 1}. In the subspace of \s) the Hamiltonian (6) 
reduces to an effective qubit Hamiltonian 



H{s) 



(8) 



with an effective external "magnetic field" b{s) 



^^^j S„s„. H{s) has eigenvalues ±f}(s) = 



±^JjU!'^ + b'^{s) with corresponding eigenstates propor- 
tional to ±[0.T 5t^]|-)- 



IV. PURE INITIAL STATE OF SPINS 

I open discussion of decoherence with an example 
where the initial state is 

1^(0)) = (a|+)+/3|-)) |li,...,M . (9) 

Each spin is initially in the +1 eigenstate of its ct^ . 
More general discussion is postponed to Section VIII, 
where I consider an ensemble of pure initial states. How- 
ever, average over the ensemble will give the same results 
as the present example so all the conclusions of this Sec- 
tion are also valid for the ensemble of spin states. 
The Hamiltonian (6) evolves the initial state into 



+ 



IP + {\uj - 

ab-l3{^cu + n) 



(10) 
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B{t, s) = -e 



ab - 



n) 



fe2 + (iw - 

ab — (3{^co 
2'^' 62 + (ia; + f2)2 



(11) 



This state is an entangled state of the qubit and the 
spins: the state of the qubit A{t,s)\+) + B{t,s)\—) de- 
pends on the state of spins \s). The state of the qubit can 
be described by reduced density matrix p{t) obtained by 
taking a partial trace over the spins: 



p{t)^T:rsm))m)\ = ^Y. 



2JV 



AB* 
BB* 



(12) 



This matrix form is given in the basis of |±) states of 
the qubit. p{t) is in general a mixed state as can be 
most conveniently measured by its "quadratic entropy" 
S{t) = 1 - Trp2(t) e [0, i] introduced in Ref. [18]. The 
initial state is pure, S{0) = 0, but in general S{t) grows 
as the interacting qubit and spins get entangled and the 
state of the qubit is loosing its initial coherence. 

The most right hand side of Eq.(12) can be formally 
interpreted as an average over all possible states of spins, 
or over all possible sets s of N independent random vari- 
ables s„. In the limit of large N the effective magnetic 
field b{s) = p + ^n=i ^nSn bccomcs a gaussian variable 
with a mean of fi and a variance of NB"^ = ^^=i ^n- 
In this limit the sum in Eq.(12) can be approximated by 
an integral. 



Pit) 



/+00 
db 
-OO 



V2ttNB^ 



A*A 
A*B 



AB* 
BB* 



(13) 



This density matrix is formally an average over different 
values of the static magnetic noise b generated by dif- 
ferent static spin configurations. The mixed state of the 
qubit results from an average over difi^erent pure states 
of the qubit, each of them evolved from the same initial 
state but with a different value of b in the Hamiltonian 
(8). However, as the derivation of Eq.(13) demonstrates, 
this is not an average over different realizations of the 
experiment with different static configurations of nuclear 
spins, but it is an average over a superposition of different 
states of spins present in a single realization of the exper- 
iment. This is genuine decoherence due to entanglement 
with the spin bath. 

In the following sections special cases of this exact solu- 
tion are worked out in more detail. I begin with the sim- 
plest case when the frequency of Rabi oscillations a; = 0. 



V. NO RABI OSCILLATIONS 

The density matrix (13) becomes particularly simple in 
the absence of flux oscillations when we can set Cl{s) = 
b{s) and 



A = a cos bt — i(3 sin bt , 
B = (3 cos bt — ia cos bt . 

The products A* A or A*B in the density matrix (13) con- 
tain oscillatory terms proportional to e^^*''* _ xhe average 
over b in Eq.(13) dephases these oscillations to zero, 



/+00 
db- 
-oo 1 



after decoherence time 
To = 



^±2int^-2NB^t'' 



2\/]VB2 ■ 



(14) 



(15) 



It is instructive to write the density matrix (13) in the 
basis of (Tx-eigenstates 



P{t) = 

|a + /3|2 
(a*+/3*)(a-/?)e" 



-2ifj,t 



, (a + /3)(a*-/?*)e 



T+2iij,t 



The off-diagonal coherences between the Uz-eigenstates 
are destroyed after the decoherence time of tq when the 
density matrix becomes diagonal in this basis. This is 
not quite surprising because the cr, cigcnstatcs arc eigen- 
states of both the interaction Hamiltonian V and the 
qubit Hamiltonian (w = 0) - they are the ideal pointer 
states of Ref. [19]. Quadratic entropy of the state (17) 
grows like 



S{t) 



1 



a + /3|2|a-/?l2 1 



(18) 



The only case when the state of the qubit remains pure, 
or S{t) = 0, is when the initial state is one of the u^- 
eigenstates {a = ±/3). This is another fundamental prop- 
erty of the pointer states. 



VI. FAST RABI OSCILLATIONS 

When Rabi oscillations are much faster than the inter- 
action with the spin bath, then we have a small parameter 



« 1 



(19) 



The matrix elements in the density matrix (13) contain 
oscillatory terms which can be approximated as 



(20) 



with a;2 = a;2 -|- Ap!^. In this expansion I use the assump- 
tion (19) that ~ e «; 1. 
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Assuming further that the bias /U is weak as compared 

to UJ, 



n"^ < 



(21) 



the density matrix can be approximated by its leading 
order term in both /i and {b — fj,): 



Pit) = 



a* (3z{t)e 



iiot 



aP*z*{t)e 



with the coherence function 



Z{t) = 



(1-^1) 



1/2 



(23) 



This function comes from an average of the approximate 
(20) over the gaussian variable b. z{t) has two time scales 



T = — < To , 

e 

V = J2^ < To 



(24) 

(25) 



both of them arc much longer tlian the dccohcrcncc time 
To in the absence of Rabi oscillations. The coherence 
function z{t) determines both the entropy growth 



S{t) = 2\am{l-\z{t)\') 



(26) 



and coherent flux oscillations 



F{t) = Tip{t)a, = -[z{t)e'''' + c.c.] . (27) 

starting from the initial +1 cigcnstatc cr^. The exact 
solution (22) can be best imdcrstood in limiting cases. 

For times i 3> t the coherence function can be ap- 
proximated by z{t ^ t) = -^/j expz ( J 
coherence decay is characterized by the power laws: 



J and the 



S(t»T) = 2|an/3|2 
F{t » t) = 



(-7) 



— cos 

t 



(28) 
(29) 



Apart from the shift in the frequency of oscillations, in 
the regime of t t the coherence decay does not depend 
on T^. This power law decay of flux oscillations was also 
derived in Ref. [17]. Both the entropy growth and the 
decay of flux oscillations are due to the decay of the same 
coherence function z{t) in the density matrix (22). Thus, 
contrary to Ref. [17], it is not possible to see oscillations 
after decay of quantum coherence: flux oscillations are 
coherent flux oscillations. 

When f <C T the coherence function \s z{t <^ t) = 

exp ^— and the decay of coherence is gaussian: 



5(t<T) = 2|a|2|/?|2(l-e 1 , (30) 



(31) 



F{t < t) = e coswi . 



When the bias /i is stronger than the influence of spins, 
NB"^ <C /x^ <C , then -^t and the gaussian decay of 
coherence is completed before crossover to the power law 
(22) decay after t. For a weak bias, /i^ <§; NB^ <C uP , the 
coherence decay is a power law. The coherence decays 
most slowly when the bias fjb = Q. The power law (29) 
and gaussian (31) decays are consistent with the decays 
derived in Ref. [20] from a model with a non-markovian 
external noise, compare also Ref. [22]. Here, these re- 
sults follow from a microscopic description of quantum 
decoherence in a bath of non-interacting spins. 

After decoherence the density matrix (22) becomes di- 
agonal in the basis of da^-eigenstates. What is more the 
entropy remains zero, S{t) = Q, only when a = or /3 = 
i.e. when the initial state of the qubit is an eigenstate of 
Gx- This is not quite surprising because when Hq dom- 
inates over V and Hg, then the eigenstates of Hq are 
expected to be the pointer states [21]. 



VII. EARLY TIME PARTIAL DECOHERENCE 

The density matrix (22) is the leading order term in 
the expansion of the exact density matrix (13) in powers 
of ^ ~ e. In this Section I include higher order terms. 
For the sake of simplicity I consider only the optimal case 
of zero bias when p, = Q. 

The products A*B etc. in the exact density matrix 
(13) are expanded up to second order in powers of ^ . The 
oscillatory factors e^*^^* are approximated as in Eq.(20). 
After these two approximations the gaussian integral over 
b in Eq.(13) gives the approximate density matrix accu- 
rate up to second order in e 



Pit) = 

with matrix elements 
P{t) 



P 



C* 



C , l-P 



\a\^- 



(32) 



(33) 



e^H^{t)+c.c: 



C{t) = a*l3z{t)e"^* + 



\p (a*/3 + c.c.) \l-e^^z^{t)\ 



,iu)t „3 



it) 



■ c.c. 



(34) 



Here P (or 1 — P) is a probability to find the qubit in the 
state 1+) (or |— )) and C is quantum coherence between 
these states. The z{t) is the function (23) with zero bias 
/U = or T„ = 00. 
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For /i = the function z{t) has only one timescale t. 
The regime t ^ t has been considered in the proceeding 
Section. In the opposite regime of c^arly time we can 
approximate z{t ^ r) = 1 in Eqs. (33,34) and get the 
entropy 

S{t<^T) = (35) 
£^(1 - cosut) {\a\^ + |/3|* - [{a*(3fe"^* + c.c] } (36) 

No matter what is the initial state, defined by a and /?, 
this entropy does not remain zero but fluctuates with an 
amplitude proportional to e^, sec Fig.l. This fluctuating 
entropy means partial loss of coherence. 

This partial early time decoherence also shows in flux 
oscillations. The initial state — P — leads to flux 
oscillations 

F{t) = Tr a^p{t) = + (1 - c^) costot . (37) 

Unlike the coherent flux oscillation between +1 and — 1 
in Eq.(2) this flux oscillates between +1 and — 1 + 2e^, 
see Fig. 2. The amplitude of the coherent oscillation is 
reduced due to the partial decoherence. The oscillation 
is also biased towards the initial positive value of F. Sim- 
ilar effect was also predicted recently in the spin-boson 
model [23] . Here it follows from the simple and explicitly 
solvable spin-spin model. The effect seems to be generic 
feature of a non-markovian environment. 

Fast Rabi oscillations make the time of full decoherence 
r much longer than the decoherence time tq without Rabi 
oscillations. When the flux oscillates many times before 
the spins can learn its orientation, then it takes the spins 
much longer time to learn the state of the qubit than 
in the static case. On a long timescale covering many 
flux oscillations the oscillating magnetic field measured 
by the spins effectively averages out to zero. This is the 
key idea of the quantum bang-bang control introduced in 
Ref. [24]. However, this argument does not apply before 
the first period of oscillation is completed. On this short 
timescale the magnetic field of the qubit is effectively 
static and during this short time the spins have a chance 
to get a rough idea about its orientation. The spins more 
or less realize that the initial flux in Eq.(37) is and 
this fuzzy knowledge biases the following oscillations in 
F(t) towards this positive value. The state of the qubit 
is partially "collapsed" towards the state with positive 
flux. This rjualitative argument becomes more substan- 
tial when we take as an initial state the equal superposi- 
ton of states with opposite flux |-|-) (or a = = 0) 
and compare the two expressions for entropy, one valid 
for fast Rabi oscillations (36) at an early time long be- 
fore the first period of oscillations is completed i.e. for 
<C 1, and the other in the absence of any oscillations 
(18). For such early times both formulas are the same, 

S{^t « 1) = , (38) 



demonstrating that the early time partial decoherence 
takes place on the timescale tq as if there were no Rabi 
oscillations at all. 

Fast Rabi oscillations postpone full decoherence after 
r but they do not prevent partial decoherence already 
after tq. This effect may affect scallability of the su- 
perconducting flux qubit technology. The partial loss of 
coherence ~ is small in a single qubit, but it scales 
together with the number of qubits A^qubits as A^qubits^^ 
and limits the number of qubits to 

A^qubits < ^ • (39) 

Obviously, this upper bound does not apply when quan- 
tum error correction is performed much faster than tq. 

VIII. ENSEMBLE OF SPINS STATES 

So far the initial state of spins was a special pure state. 
In this Section I assume the unbiased initial ensemble of 
all possible spins' states described the unit density matrix 
of spins Is- The initial state of the qubit and the spins 
is 

PQs(O) = (a|+) + /3|-))(a*(+| +/?*(-[) ® Is- 

(40) 

The unit matrix Is can be expressed in any basis of spins 
states but a representation Is = -^J^sl^i^ most di- 
rectly leads to the reduced density matrix of the qubit 

p{t) = Tr sPQs{t) = ^ X] ( A*B BB* ) ' ^^^^ 

This density matrix is identical with the density matrix 
(12) obtained from the initial pure state in Eq.(9). All 
the results derived from the pure state pass without any 
modification to the ensemble. The initial state I5 can be 
also interpreted as a high temperature thermal state of 
nuclear spins. 

IX. ESTIMATES 

In this Section I attempt to make some numerical es- 
timates. For example, in the qubit design considered 
in Ref. [16] the parameters are estimated as: N = 10*, 
B = 10^' '^s"^ and co = lO^s^^. These numbers lead to 
TO = 10-8--9s, = IO-2 -O, and r = IQ-^ -Ss. In 
the best case of = 1G~^ one can see u)t = = 10^ 
of coherent fiux oscillations before the decoherence time 
T. Thus it is possible to see hundreds of oscillations as 
in the experiments [7]. For the same parameters the 
fast partial decoherence limits the number of qubits to 
-^qubits ^ 10^- These numbers leave, however, an en- 
couraging room for current experiments. 
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X. CONCLUSION 

The environment of spins is a witness to the quantum 
state of the superconducting flux qubit. The qubit gets 
entangled with spins after a decoherence time tq. Co- 
herent flux osciUations driven by the qubit Hamiltonian 
postpone full decoherence after a much longer time r, but 
they do not prevent partial decoherence already after tq. 
This partial decoherence biases the flux oscillations to- 
wards the initial orientation of the flux. 
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FIG. 1. Entropy S{t) according to the exact Eq.(13) for 
w = 1, /i = 0, e = 0.1, and a = P = The timescale for full 
decoherence is r = 400 and the early time partial decoherence 
sets in on the timescale tq = 10. The inset shows the entropy 
at early times. The non-vanishing early time entropy means 
partial decoherence. 
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FIG. 2. Flux F{t) according to the exact Eq.(13) for w = 1, 
/i = 0, e = 0.2, and — /3 = As a result of the fast partial 
decoherence the flux oscillation is biased towards the initial 
F{0) = +1. 
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